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Abstract 

We consider the weak decays of a B meson to final states that contain a S-wave ra- 
dially excited vector meson. We consider vector-pseudoscalar final states and calculate 
ratios of the type B — > p'n/B — > pir, B — » lo'tt/B — * uir and B — > (fi'ir/B — * cpir where 
p', to' and (f)' are higher p, uj and (ft S-wave radial excitations. We find such decays to 
have larger or similar branching ratios compared to decays where the final state p, uj 
and (ft are in the ground state. We also study the effect of radial mixing in the vector 
system generated from hyperfine interaction and the annihilation term. 
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The new data accumulating from B factories and other accelerators will include transi- 
tions to many new final states which have not been previously studied in detail; e.g. radially 
excited meson states. Many decays involve a transition from a low momentum spectator 
quark to a high momentum relativistic meson. The form factors for such transitions are ex- 
pected to be sensitive to the high momentum components of the final meson wave function, 
and therefore to favor radially excited states. The data on B decays to these states will thus 
provide important new information, particularly for the form factors to the radially excited 
states and probe the high-momentum tails of their wave functions. 

In this paper we calculate predictions for the ratios 

R p+ = BR(B° -f 7i-p + ')/BR(B° -> 7r~p + ) (1) 

R p0 = BR(B~ -> 7i-p°')/BR(B- -> Tr^p ) (2) 

R u = BR(B~ -> ir-u')/BR{B- -> tTw) (3) 

R+ = BR(B s ^-k q <P')/BR(B s ^h q <P) (4) 

where p',uj' and 0' are the radially excited states. Most studies of two body nonleptonic 
B decays have concentrated on processes of the type B — ► M\Mi where both M\ and M2 
are mesons in the ground state configuration. Nonleptonic decays, where one of the mesons 
in the final state containing the spectator quark is a radially excited state, are expected to 
have larger or similar branching ratios compared to decays where the final state contains the 
same meson in the ground state. This is easily seen in a simple model in which B — > ttM 
and M is a simple flavor eigenstate with no flavor mixing beyond isospin. We follow the 
inactive spectator approach used [|TJ to treat B decays to charmonium in which the spectator 
quark does not participate in a flavor- changing interaction and later combines with a light 
antiquark to make the final light meson as shown in Fig. [I]. The decay amplitudes are 
then described as the product of a b-quark decay amplitude and a hadronization function 
describing the combination of a quark- antiquark pair to make the final meson. Neglecting 
the relative Fermi momentum of the b quark and the spectator quark in the B meson, the 
quark transition for the processes in Eq. (1-3) is 

b -> Tr~(p)u(-p) (5) 
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where the b quark is at rest and p denotes the final momentum of the n 
in Eq.[|, the quark transition is essentially similar to the one above 

b — > 7r (p)s(—p) 

where now p denotes the final momentum of the 7r°. 



B(0) 




q(0) 

Figure 1: Factorization for the decay B — > Mir. 



'. For the process 
(6) 



Concentrating on the processes in Eq. (1-3) the transition matrix for the full decay has 
the form|l| 

n-(p)M(-p)\ T \B) = (M(—p) \ F \u(-p)q(0)) ' {^{p)u{-p))\ W \b) (7) 

where T denotes the transition matrix for the hadronic decay which factors, as shown in 
Fig. [|, into a weak matrix element at the quark level denoted by W and a recombination 
matrix element denoted by F. This latter matrix element describes the transition of a quark 
with momentum — p and an antiquark with zero momentum to make a meson with momen- 
tum —p. Radial excitations could be favoured over ground states, if the final momentum p 
is large, since the radial excitations are expected to have higher kinetic energies. 

An alternative but equivalent way of understanding this effect is to note that for a heavy 
b quark one can write @ 



7T-(p)M(-p) T\B) = {M{-p)\J 1 ,\B)-U-{p) J£|0> 



(8) 

where are currents that occur in W = J\ x J 2 . The transition matrix element for 
the hadronic decay can then be written in terms of B — > M form factors and the pion 
decay constant. The form factors can be expressed as overlap integrals of the B and the 
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M meson wavefunctions. When M is a light meson, with a mass much smaller than the B 
meson, the overlap integrals get contributions mainly from the high momentum components 
of the meson wavefunctions. It is now clear that for a radially excited meson M', which has 
more higher momentum components, the overlap integrals will be enhanced compared to the 
overlap integral with a ground state meson M. Consequently the B — > M' form factors are 
likely to be enhanced relative to the B — > M form factors which would then translate into 
higher branching ratio for B — > M'tt relative to B — > Mir. 

Our discussion above assumed the physical states to be pure radial excitations. However, 
additional interactions can mix the various radial excited components. For instance hyperfine 
interactions can mix radial excitations with the same flavor structure and so in general, in 
the p, to and <fi systems, the various physical states will be admixtures of radial excitations 
P. f|. Flavor mixing in the vector system is known to be small but is important in the 
pseudoscalar sector. We will consider the pseudoscalar case in a different publication ||. To 
make quantitative predictions we use constituent quark wave functions with several potentials 
to test the dependence of the results on the confining potential. We shall see that the effects 
of the potential dependence and mixing are small so that the results are reasonably robust 
and are not seriously dependent on the fine details of the model. 

Even though our discussion has so far only included vector-pseudoscalar final states we 
can also consider vector- vector final states such as B — > pp or B — > J/ipK*. However vector- 
vector final states are complicated since different partial waves are present. Our purpose 
here is to demonstrate the effects of radial mixing in the simpler physical system of the 
vector-pseudoscalar final state. If the effects of radial enhancements are observed in the 
vector-pseudoscalar case we would expect them to be also present in the vector-vector final 
state. 

We will first review the study the masses and mixing in the vector meson sector following 
the simple nonrelativistic approach in Ref@, f|]. To obtain the eigenstates and eigenvalues 
in the vector meson system we diagonalize the mass matrix 

< q' a q' b , n'\M\q a q b , n> = 5 aa '8 w 8 nnl {m a + m b + E n ) + 5 aal 5 w s a ■ s b i[> 

m a m b 

where s a ^ and m a ,b are the quark spin operators and masses. Here n = 0, 1, 2 and the basis 
states for the isovector mesons are chosen as \N, I — 1, 1 3 — 1} — — uctj, \N, I — 1, 1% — 0) = 
uu — dctj /\/2 and \N,I = 1, 1$ = —1) = \du) where 1,1^ stand for the isospin quantum 
numbers. In the above equation E n is the excitation energy of the n th radially excited state 
and B is the strength of the hyperfine interaction. 
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To begin with, we use the same harmonic confining potential as well as the other param- 
eters used in Ref|4j] to obtain the eigenstates and eigenvalues for the mass matrix in Eqn. || 
The various parameters used in the calculation are the constituent masses, m u = = 
0.350 GeV, m s = 0.503 GeV, the angular frequency, u =0.365 GeV and b = B/m 2 u =0.09. 
The eigenvalues and eigenstates for the p system with a harmonic potential are shown in in 
Table. [l|. To see how this result changes with a different confining potential we use a power 

Table 1: Eigenvalues and Eigenstates for the p system-Harmonic potential 



Harmonic 


iV 


Ni 


N 2 


p(0.768) 


0.990 


0.124 


-0.066 


p(1.545) 


0.108 


-0.973 


0.204 


p(2.370) 


0.089 


-0.195 


0.977 



law potential V(r) = Xr n . We will use a linear and a quartic confining potential and 
compare the spectrum with that obtained with a harmonic oscillator potential. To fix the 
coefficient A we require that the energy eigenvalues of the Schrodinger equation are similar 
in the least square sense with the energy eigenvalues used in RefQ. So for example, for the 
linear potential, we demand that 

F = ^2,{E n {harmonic) — E n (linear)) 2 

n 

is a minimum. This fixes the constant A in V(r) = Xr and we obtain the eigenvalues and 
eigenstates in Table. ||. We follow the same procedure for the quartic potential and obtain 

Table 2: Eigenvalues and Eigenstates for the p system-Linear potential 



Linear 


N 


N x 


iV 2 


p(0.775) 


0.992 


0.112 


-0.053 


p(1.515) 


0.104 


-0.986 


0.130 


p(2.260) 


0.066 


-0.122 


0.990 



the eigenvalues and eigenstates in Table. ||] We observe from Tables. (1-3) that the mass 
eigenstates and eigenvalues of the p system are not very sensitive to the confining potential 
and the radial mixing effects are small. 



5 



Table 3: Eigenvalues and Eigenstates for the p system- Quartic potential 



Quart ic 


iV 


Nt 


N 2 


p(0.759) 


0.988 


0.129 


-0.077 


p(1.567) 


0.103 


-0.955 


0.278 


p(2.370) 


0.11 


-0.267 


0.957 



To obtain the eigenstates and eigenvalues in the u — system we diagonalize the mass 
matrix 

< q' a c? b , n'\M\q a q bl n> = 5 aa >8 w 5 nn/ (m a + m b + E n ) + 5 aa >5 w s a ■ s b ^ 

m a m b 

+ Sab6a'b>^—MWn'(0) (10) 

m a m b 

This has a similar structure as the p system but now we have the additional annihilation 
interaction with strength A that causes flavor mixing. |3|, [7]. 

Diagonalizing the mass matrix in Eqn. 1C , with the basis states \N) = uu + dctj / \[2 and 
\S) = \ss), we obtain the eigenvalues and the eigenstates of the uj — <fi system. We use the 
same value for the hyperfine interaction as used for the p system. For the linear potential 
we obtain with B = 0.09m^ and A = 0.005m^ the eigenstates and eigenvalues in Table. |]. 
For the harmonic potential we obtain with B = 0.09m^ and A = 0.015m^ the eigenstates 



Table 4: Eigenvalues and Eigenstates for the uj — 4> system- Linear potential 



Linear 


iV 


N x 


N 2 


So 


Si 


s 2 


cj(0.782) 


0.991 


0.123 


-0.058 


-0.014 


0.004 


-0.002 


0(1.05) 


0.012 


0.011 


-0.004 


0.997 


0.071 


-0.034 


w(1.52) 


-0.113 


0.982 


0.144 


-0.006 


-0.034 


0.004 


0(1.66) 


0.007 


-0.030 


-0.014 


0.068 


-0.994 


-0.077 



and eigenvalues in Table. [|. For the quartic potential we obtain with B = 0.09m^ and 
A = 0.023m^ the eigenstates and eigenvalues in Table. || As in the p system we find the 
mixing to be insensitive to the confining potential and the effects of radial mixing to be 
small. We also find, as expected, a small value for the annihilation term in the fits to the 
masses. 
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Table 5: Eigenvalues and Eigenstates for the u — system- Harmonic potential 



Harmonic 


iV 


iVi 


N 2 


So 


Si 


s 2 


u;(0.783) 


0.984 


0.154 


-0.081 


-0.033 


0.011 


-0.007 


0(1.05) 


0.026 


0.029 


-0.011 


0.994 


0.089 


-0.048 


w(1.57) 


-0.126 


0.948 


0.256 


-0.008 


-0.139 


0.010 


0(1.68) 


0.025 


-0.12 


-0.07 


0.082 


-0.976 


-0.143 



Table 6: Eigenvalues and Eigenstates for the u — system- Quartic potential 



Quart ic 


N 


Ni 


N 2 


So 


Si 


s 2 


cj(0.783) 


0.980 


0.163 


-0.096 


-0.049 


0.012 


-0.009 


0(1.05) 


0.041 


0.034 


-0.015 


0.991 


0.100 


-0.060 


w(1.58) 


-0.122 


0.932 


0.322 


-0.010 


-0.113 


0.006 


0(1.7) 


0.022 


-0.089 


-0.067 


0.086 


-0.968 


-0.207 



We now use these wavefunctions to predict the ratios in Eq. (1-4). These decays are 
dominated by diagrams which satisfy the inactive spectator approach and are treated 
with Eq.(7-8). Some of the diagrams which violate this assumption; e.g. penguin and 
annihilation contributions, may not be as negligible here as in the charmonium case treated 
in ReffllJ for the decays to the ground state configurations. But they are expected to have 
much smaller form factors for radial excitations. Therefore it is reasonable to neglect them 
for this preliminary investigation of the order of magnitude of these ratios. Note that for 
the B s -> 7T°0 decay the QCD penguin is isospin forbidden, the annihilation contribution is 
OZI forbidden and the electroweak penguin is also described by Eq. [7|. 

We obtain, using factorization for the nonleptonic amplitude, 



(p +, \ 


U7"(l 


-7s)b 


B°) 


) (tt 1 djn(l 


-Tb)«|0> 


(p + \ 


M7^(l 


-lh)b 


B°) 


(71- 1 (?7 M (1 


- l5 )u\0) 
















p3 

P 











(11) 



where P is the magnitude of the three momentum of the final states and the form factor A' 
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is defined through 



(VflA^Pj = (M l + M f )A 1 



+ 2M f A Q e -^q„ 



e .q 
q 2 



A, 



e*.q 



'Mi + M 



f L 



(Pi + Pf)n - 



M? -M 



'12) 



where A^ is the axial vector current. Similarly we obtain 



R p o 



(p°'\ 


U7^(l 


-75)6 


BP) (tt- 


1 ^(1 


-75)«|0) 


(p°\ 




-75)6 


B°\ (tt- 


^7/i(l 


-75^10) 




2 p3 























(13) 



Finally, 



(u/| 


M7"(l 


-75)6 


5°) 




- 7 5 )m |0) 




■U7^(l 


-75)6 




(7r-|d7 M (l 


- 7> |0) 


A>« 


2 p3 
V 
p3 











(14) 



s 7 ^(l -75)6 


B a ) (vr 


«7m(1 -7s)« |0> 


(0| S 7 A*(1 - J 5 )b 


5°) (tt° 


«7m(1 - 7s)^ |0) 



An 



An 



p3 

Pi 



(15) 



To calculate the ratios we need the form factor Aq. Note that the wavefunctions for the 
various vector meson states are not enough to calculate non leptonic decay amplitudes. 
In particular, with the factorization assumption for non leptonic decays, the calculations 
of decay amplitudes require the matrix elements of current operators between the physical 
states. These matrix elements can be expressed in terms of form factors and decay constants. 
In this work we use a constituent quark model(CQM) model for the form factors || that 
incorporates some relativistic features and is relatively simple to adapt to the case of radially 
excited states. In this model the form factor A is given by 



An 



Z[h- 



M_ 
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where 



Z 



\M,M r .\I. 
Ml - q 2 

r3„ 



h = m s / d 3 p(f)* f (p + a)fa(p)[ 



[io? rrif 



M± = M f ±Mi 

m s q 



2m s (3 = 2 



M+ 



q + Ml-q 2 
mi + m/ 

and 0/ and fa represent the momentum space wave functions while (3 is the velocity of the 
mesons in the equal velocity frame ( also called the Breit frame or the brick wall frame) 
and rriij are the non spectator quark masses of the initial and final meson. The equal 
velocity frame in a convenient frame to calculate the Lorentz invariant form factors where 
the velocities, and f3f of the mesons with masses Mi and Mf are equal in magnitude but 
opposite in direction. We use the momentum wavefunction <pf obtained from spectroscopy 
in section. (2) while for fa we use the wave function 

fa = <P B = N B e- p2 ^ (17) 

where pp is the Fermi momentum of the B meson. In our calculations we will take Pf — 
300MeV. Note that in the analysis presented in the introduction we have neglected the Fermi 
momentum of the b quark, since pF/rrib is small. 

For transitions to higher resonant states, we use the same quark masses as those used 
in the transition of the B meson to the lowest resonant state. This is reasonable, as the 
spectator quark still comes from the B meson and therefore has the same value for its mass 
irrespective of whether the final state is in the lowest or the first excited state. The values 
for the masses of the non spectator masses are taken to be the same as those used for 
spectroscopy. However for the calculation of the velocity f3 and hence a defined in Eqn. [16| 
we use the physical mass of the higher resonant state. In Table. [7| we give our predictions for 
the various ratios defined above. We find that the transitions to higher excited states can 
be comparable or enhanced relative to the transitions to the ground state. From Table. 
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Table 7: Ratios of branching ratios for different confining potentials 



Ratio 


Linear 


Quadratic 


Quart ic 


R p + 


2.3 


2.0 


1.9 


R p o 


2.3 


2.0 


1.9 


Ruj 


3.5 


2.5 


1.7 


R<p 


6.7 


6.2 


5.2 



we see that the ratios of branching ratios are slightly sensitive to the confining potential and 
the ratios of branching ratios increase as we go from the quartic to the linear potential. This 
is because the wavefunction for the linear potential has a longer tail and hence more high 
momentum components than the wavefunction for the quadratic and the quartic potentials. 
The wavefunction for the quadratic potential, has in turn, a longer tail and hence more 
high momentum components than the wavefunction for the quartic potential. Hence we 
would expect the hierarchy {A ) Hnear > (A ) quadraUc > {A ) quartic and a similar one for the 
radially excited states (A' ) Unear > (A' ) quadratic > {A' ) quartic where A and A' are the form 
factors for the transition of B to the ground state and the first radially excited state of the 
meson M. We see from from Table. [5] that this hierarchy is maintained for the ratios of 
form factors and so we have (A' Q / 'A Q ) Unear > (A' Q / A ) quadratic > (A' Q / A ) quartic . Note that the 
ratio of form factors also depend on the choice of the Fermi momentum of the B meson, as 
a smaller(larger) Fermi momentum would make the form factors more(less) sensitive to the 
tail of the wavefunction of M, as well as mixing effects in the wavefunction of the meson 
M. The effect of mixing between the various radially excited states and the ground state is 
generally small. 

We observe in Table. [7| that there can be a large enhancement for R^. This decay is 
suppressed in the standard model. One can get a rough estimate of the branching ratio for 
B s — > (f)TT° using factorization as 



BR[B S -> 07T°] _ 1 
BR[B° -»• P+7T-] ~ 2 



V ub V u * s ( Cl + c 2 /N c ) - V tb V t * s 3c 9 /2 



0.02 



V ub V: d (c 2 + Cl /AQ 

where we have neglected form factor and phase space differences between B s — > <f)n and 
B° — > p + it~. The Wilson coefficients q can be found in Ref|8| while V u b,V us ,V u d,V t b and 
V ts are the various CKM elements UI0" |. Using the measured BR[B — ► p + 7i~] ~ 28 x 10~ 6 
]TT| we get BR[B S — > 07r°] ~ 6 x 10~ 7 . Hence the large enhancement for indicates that 
BR[B S -> 0'tt ] can be around ~4x 10" 6 . 
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Note that in the p{u) system there are two resonances, p(1450)[u;(1420)] and p(1700)[o;[1650], 
which can be identified a S-wave radial excitation(2S) and a D wave orbital excitation in 
the quark model. However recent studies of the decays of these resonances show that it 
is possible that these states are mixtures of qq and hybrid states Ref|IJ. Hence the state 
p(1450)[cg>(1420)] is interpreted as a 2S state with a small mixture of a hybrid state. We 
do not take into account such possible mixing with a hybrid state in our calculation and 
the meson masses for these excited states used in our calculation are the ones we predict in 
section 2. For the system there is only state at 0(1680) which we interpret as a 2S state 
in the absence of mixing effects. 

In conclusion we have considered the weak decays of a B meson to final states that are 
mixtures of S-wave radially excited components. We calculated nonleptonic decays of the 
type B — > p'tt/B — > pn, B — > u'it/B — ► uj-k and B — > <p'ir/B — ► 07r where p', uj' and 0' are 
higher p, uo and resonances. We found that the transitions to the excited states can be 
comparable or enhanced relative to transitions to the ground state. It would, therefore, be 
extremely interesting to test these predictions. We also studied the effect of radial mixing 
in the vector system generated from hyperfine interaction and the annihilation term; these 
turn out to be generally small. 
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